Speckle Imaging of Spin Fluctuations in a Strongly Interacting Fermi Gas 
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Spin fluctuations and density fluctuations are studied for a two- component gas of strongly inter- 
acting fermions along the BEC-BCS crossover. This is done by in-situ imaging of dispersive speckle 
patterns. Compressibility and magnetic susceptibility are determined from the measured fluctua- 
tions. This new sensitive method easily resolves a tenfold suppression of spin fluctuations below 
shot noise due to pairing, and can be applied to novel magnetic phases in optical lattices. 
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One frontier in the field of ultracold atoms is the real- 
ization of quantum systems with strong interactions and 
strong correlations. Many properties of strongly corre- 
lated systems cannot be deduced from mean density dis- 
tributions. This has drawn interest toward novel ways 
of probing cold atoms, e.g. via RF spectroscopy [ll|2], 
Bragg and Raman scattering [3 , interferometric methods 
^ 5 and by recording density correlations [6-8 . Fur- 
ther insight into quantum systems is obtained by looking 
not only at expectation values, but also at fluctuations. 
Several recent studies looked at density fluctuations, ei- 
ther of bosons around the superfluid-to-Mott insulator 
transition [9UTT]. or of a gas of non-interacting fermions 

mils]. 

In this paper, we extend the study of fluctuations of 
ultracold gases in several ways. First, we apply it to a 
two-component Fermi gas across the BEC-BCS crossover. 
Second, we implement a very sensitive way to measure 
fluctuations in the magnetization, i.e. the difference of 
the densities in the two different states. Third, we in- 
troduce the technique of speckle imaging as a simple and 
highly sensitive method to characterize fluctuations. 

Our work is motivated by the prospect of realizing 
wide classes of spin Hamiltonians using a two-component 
gas of ultracold atoms in an optical lattice [Ml [15] . An 
important thermodynamic quantity to characterize two- 
component systems is the spin susceptibility, which pro- 
vides a clear signature of phase transitions or crossovers 
involving the onset of pairing or magnetic order [16| Il7j. 
At a ferromagnetic phase transition the susceptibility di- 
verges, whereas in a transition to a paired or antiferro- 
magnetic phase the susceptibility becomes exponentially 
small in the ratio of the pair binding energy (or antifer- 
romagnetic gap) to the temperature. The fluctuation- 
dissipation theorem relates response functions to fluctu- 
ations, consequently the spin susceptibility can be de- 
termined by measuring the fluctuations in the relative 
density of the two spin components. 

In our experiment spin fluctuations create phase shifts 
of (detuned) imaging light that vary randomly in space; 
we measure them by imaging the resulting speckle pat- 
terns. When imaging atom clouds, one usually relates the 
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FIG. 1: Simulation of propagation effects after light has 
passed through a Poissonian phase noise object. Shown are 
the variance measured in the amplitude or in-phase quadra- 
ture (black line) and the out-of-phase quadrature (gray line) 
as a function of defocus distance, for an imaging system with 
a numerical aperture of 0.14. Within a distance less than 5 
percent of our cloud size, noise becomes equally distributed 
between the two quadratures and the variances in transmis- 
sion and phase- contrast images become the same. (Top inset) 
For small phase fluctuations, an in-focus phase noise object 
gives no amplitude contrast, but when it is out of focus it 
does. (Bottom inset) Sample intensity patterns for a defo- 
cused phase object. 



transmitted light intensity with absorption and its phase 
with dispersion fT8|. This is different in our method. Spin 
and density fluctuations occur on all spatial scales down 
to the interatomic separation, and their observation is 
limited by the maximum resolution of the imaging sys- 
tem. For typical atom clouds in 3D (in contrast to 2D ex- 
periments [OHTTj or very small clouds in 3D [12]) the cloud 
size is larger than the Rayleigh range corresponding to 
the maximum resolution of the imaging systems. There- 
fore, for the smallest resolvable fluctuations, the depth of 
field is smaller than the sample size and, necessarily, the 
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recorded image is modified by propagation effects. Prop- 
agation mixes up amplitude and phase. This can be eas- 
ily seen in the case of a phase grating, which creates an in- 
terference pattern further downstream. Conversely, due 
to Fermat's principle, for an amplitude object rephasing 
takes place only in conjugate planes, and therefore out of 
focus there is a phase-contrast signal. Similar physics is 
responsible for laser speckle when a rough surface scat- 
ters light with random phases [19], and occurs when a 
Bose-Einstein condensate with phase fluctuations devel- 
ops density fluctuations during ballistic expansion [20 , or 
when a phase-contrast signal is turned into an amplitude 
signal by deliberate defocusing [21]. 

In our experiments, we use off-resonant imaging light 
due to the high optical density of the cloud. Absorption 
decreases inversely proportional to the detuning squared, 
whereas the dispersive signal falls off more slowly, in- 
versely proportional to the detuning. Therefore, the 
dominant signal is initially dispersive, but is converted 
to an amplitude signal during propagation. Simulations 
using Poissonian noise confirm this picture (Fig. [T]): Af- 
ter a short propagation distance, the power of phase and 
amplitude noise are equal, independently of whether the 
noise was purely dispersive or absorptive before propaga- 
tion. This feature of speckle makes our imaging technique 
both simple and robust. It is insensitive against defocus- 
ing, and allows us to image fluctuations of the real part of 
the refractive index (i.e. a phase signal) without a phase 
plate or other Fourier optics. 

The experiments were performed with typically 10^ ^Li 
atoms in each of the two lowest hyperfine states |1) and 
1 2) confined in an optical dipole trap oriented at 45° to 
the imaging axis with radial and axial trap frequencies 
ujr = 27ix 108.9(6) s-i and cj^ = 27r x 7.75(3) s-^ For 
the samples imaged at 52 7G, the sample preparation was 
similar to that described in [13], with a temperature of 
0.14(l)Ti?. The samples imaged at other magnetic fields 
were prepared in a similar fashion, except that evapo- 
ration was performed at lOOOG to a final temperature 
of T = 0.13(l)Tp before ramping the magnetic field over 
1.5s to its final value. The temperature at lOOOG was de- 
termined by fitting a noninteracting Thomas-Fermi dis- 
tribution in time of flight. The temperatures at other 
points in the crossover were related to that value as- 
suming an isentropic ramp, using calculations presented 
in [2^. Using this method we obtain temperatures of 
0.13(1)Tf at 915G, 0.19(1)Tf at 830G, and 0.19(3)Tf 
at 790G where additional evaporation was performed to 
achieve a central optical density similar to that at the 
other magnetic fields. The extent of the cloud along 
the imaging direction was 135/im, much larger than the 
Rayleigh range of 8/im for our imaging system with a 
numerical aperture of 0.14. 

The superfluid to normal phase boundary was deter- 
mined by measuring condensate fraction (Fig. [4| using 
the standard magnetic field sweep technique [2^ [24] . 



For this, the magnetic field was rapidly switched to 
570G to transfer atom pairs to more deeply bound pairs 
(molecules) which survive ballistic expansion. For reso- 
nant imaging of the molecules, the field was ramped back 
to 790G over 10 ms. The condensate fraction was deter- 
mined by fitting the one- dimensional density profiles with 
a bimodal distribution. 
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FIG. 2: (Color online) (Top) Example speckle noise im- 
age, with white box indicating analysis region. (Bottom) 
Noise data for noninteracting (left) and resonantly interacting 
(right) cold clouds, showing A?_ (black dots) and (grey 
dots). Solid lines are Gaussian fits to the data, and dotted 
lines are expected full Poissonian noise for the correspond- 
ing quantities based on density determined from off-resonant 
absorption. 

As previously described, propagation converts spatial 
fluctuations in the refractive index into amplitude fluctu- 
ations on the detector. For different choices of the probe 
light frequency, the two atomic spin states will have dif- 
ferent real polarizabilities and the local refractive index 
will be a different linear combination of the (line-of-sight 
integrated) column densities rti and 712. To measure the 
susceptibility we choose a probe light frequency exactly 
between the resonances for states |1) and |2), so that the 
real polarizabilities are opposite and the refractive index 
is proportional to the magnetization (ni — 712). The in- 
tensity fluctuations on the detector after propagation are 
consequently proportional to the fluctuations in magneti- 
zation. Since a refractive index proportional to (ni +77,2) 
occurs only in the limit of infinite detuning, we measure 
the fluctuations in the total density by exploiting the fact 
that the fluctuations in total density can be inferred from 
the fluctuations in two different linear combinations of 
ni and 77,2. For convenience, we obtain the second linear 
combination using a detuning that has the same value, 
but opposite sign for state |2), and therefore three times 
the value for state |1). With this detuning, we record 
images of the fluctuations in (77,1 /3 + 77,2). 

In principle this information can be obtained by tak- 
ing separate absorption images on resonance for states 
|1) and |2). However, the images would have to be taken 
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on a timescale much faster than that of atomic motion, 
and for there would be increased technical noise from 
the subtraction of large numbers. The use of dispersive 
imaging has the additional advantage over absorption in 
that the number of scattered photons in the forward di- 
rection is enhanced by superradiance. As a result, for the 
same amount of heating, a larger number of signal pho- 
tons can be collected [18 . This is crucial for measuring 
atomic noise, which requires the collection of several sig- 
nal photons per atom. The choice of detuning between 
the transitions of the two states has the important fea- 
ture that the index of refraction for an equal mixture 
fluctuates around zero, avoiding any lensing and other 
distortions of the probe beam. This is not the case for 
other choices of detuning, and indeed, we observe some 
excess noise in those images (see below). At the detun- 
ings chosen, 10 percent residual attenuation is observed, 
some due to off-resonant absorption, some due to disper- 
sive scattering of light out of the imaging system by small 
scale density fluctuations. The contribution to the vari- 
ance of the absorption signal relative to the dispersive 
signal scales as (2r)^/^^ ^ 0.006 and can be neglected in 
the interpretation of the data. 

The noise analysis procedure was nearly identical to 
that performed in [13 . A high-pass filter with a cutoff 
wavelength of 13 /im was applied to each image of the 
cloud to minimize the effect of fluctuations in total atom 
number. Then, for each pixel position, the variance of the 
optical densities at that position in the different images 
was computed. After the subtraction of the contribution 
of photon shot noise, the resulting variance image reflects 
the noise contribution from the atoms. 
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FIG. 3: (a) The ratio x/^? (b) the normalized susceptibil- 
ity x/xo, and (c) the normalized compressibility k/kq in the 
BEC-BCS crossover. The variances derived from sequences 
of images are converted into thermodynamic variables using 
the measured temperatures and a calibration factor deter- 
mined from the noninteracting gas. The vertical line indicates 
the onset region of superfluidity, as determined via conden- 
sate fraction measurements. The curves show theoretical zero 
temperature estimates based on 1st (dotted) and 2nd order 
(solid) perturbative formulas obtained from Landau's Fermi- 
liquid theory integrated along the line of sight, and results 
from a Monte Carlo calculation (dashed) for the compress- 
ibility in a homogeneous system [25] . 



The goal of our noise measurements is to determine 
at various interaction strengths the normalized suscep- 
tibility X = x/Xo and compressibility ^ = n/n^^ where 
Xo = 3n/2£^i? and = 3/2nEF are the susceptibility 
and compressibility of a zero-temperature non interact- 
ing Fermi gas of the same total density n and Fermi en- 
ergy Ep. Before studying spin fluctuations through the 
BEC-BCS crossover, we therefore calibrate our measure- 
ment by measuring the spin fluctuations in a noninter- 
acting mixture, realized at 527G where the scattering 
length between the two states vanishes. Fig. |2] shows 
raw profiles of the difference and sum variances of the 
measured optical densities A?_ = (cA(A^i — N2))'^ and 
= {c' A{Ni/3 + A^2))^. In these relations c and c' 
are the conversion factors between number fluctuations 
and fluctuations in optical density in the specific probe 
volume V. Without interactions, Ni and N2 are un- 
correlated, and one predicts (A(7Vi - 7V2))^/(A(7Vi/3 + 
TVs))^ = 2/(1 + (1/3)2) = 1.8. The observed ratio of 
A?_/A^ = 1.56(14) reflects excess noise contributing to 
A^ due to residual systematic dispersive effects and is 
accounted for by setting c'/c — y^I.8/1.56. For high 
temperatures, the atomic noise of the non-interacting gas 
approaches shot noise; for lower temperatures we observe 
a reduction in noise due to Pauli blocking as in our pre- 
vious work f[3]. With our new method, we easily discern 
spin fluctuations with a variance of less than 10 percent 
of atom shot noise. 

The fluctuation dissipation theorem connects the vari- 
ances (A(7Vi - N2)f and (A(7Vi + TVs))^ to the suscep- 
tibility X and the compressibility H via (A(A^i — A^2))^ = 
3Ar/2 (TITf) X and (A(7Vi + N^))" = 3N/2 {T/Tp) n 
with N = Ni ^ N2 and T/Tp being the tempera- 
ture measured in units of the Fermi temperature Tp. 
Recomposing the variances from the two experimen- 
tally accessible linear combinations these relations be- 
come /S?_/Nc^ = 3/2 {T/Tf) X and 9/4 Al/Nc'^ - 
1/4 /S?_/Nc^ = 3/2 {T/Tf) n and allow us to deter- 
mine c and c' by using the 52 7G noise measurements for 
which X = ^ = 1 + ^(T/Tp)^y This analysis ignores 
corrections due to line-of-sight integration. 

Fig. [3] shows the spin susceptibility, the compressibil- 
ity, and the ratio between the two quantities for the in- 
teracting mixtures as the interaction strength is varied 
through the BEC-BCS crossover. The susceptibility and 
compressibility reproduce the expected qualitative be- 
havior: for the sample at unitarity, where the transition 
temperature is sufficiently high that a sizable portion of 
the sample is superfluid, and for the sample on the BEG 
side, where even the normal-state atoms form molecules, 
the spin susceptibility is strongly suppressed relative to 
the compressibility. This reflects the fact that the atoms 
form bound molecules or generalized Cooper pairs; the 
spin susceptibility should be exponentially small in the 
binding energy, while the enhanced compressibility re- 
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fleets the bosonie eharacter of the molecular condensate. 
At 915G and lOOOG, where the sample is above the su- 
perfluid critical temperature, the susceptibility is larger 
but still below its value for the noninteracting gas, re- 
flecting the persistence of pair correlations even in the 
normal phase of the gas. 
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FIG. 4: Measured condensate fraction as a function of di- 
mensionless interaction strength l/(kFCi)- Insets show typical 
images from which the condensate fraction was extracted by 
fitting a bimodal distribution. The dashed line is a sigmoidal 
fit to guide the eye. 



Above the Feshbach resonance, for attractive interac- 
tions, we compare our results to flrst and second order 
perturbation theory in the small parameter kpa. This ig- 
nores the instability to the superfluid BCS state at expo- 
nentially small temperatures. The perturbation theory is 
often formulated for the Landau parameters for a Fermi 
liquid [ini |26] . The susceptibility and compressibility are 
given by xo/x = (1 + F^)m/m'', kq/k. = (1 + F^)m/m*, 
where m* = m(l + is the effective mass, and Ff ^ 

Fj^ are the l-Xh angular momentum symmetric and an- 
tisymmetric Landau parameters, respectively. Although 
the experimental data are taken for relatively strong in- 
teractions outside the range of validity for a perturba- 
tive description, the predictions still capture the trends 
observed in the normal phase above the Feshbach res- 
onance. This shows that more accurate measurements 
of the susceptibility, and a careful study of its tempera- 
ture dependence, are required to reveal the presence of a 
possible pseudogap phase. 

In our analysis we have neglected quantum fluctuations 
which are present even at zero temperature [I6l[27]. They 
are related to the large-g static structure factor S{q) mea- 
sured in [22 and proportional to the surface of the probe 
volume, scaling with A^^/^ log (A^). For fluctuations of 
the total density, their relative contribution is roughly 
A/'~^/^/(T/Tp), and at most 40 percent for our experi- 
mental parameters. Attractive interactions and pairing 
suppress both the thermal and quantum spin fluctua- 
tions, but it is not known at what temperature quantum 
fluctuations become essential. 

Spin susceptibilities can also be obtained from the 
equation of state which can be determined by analyz- 



ing the average density proflles of imbalanced mixtures 
[29^. Our method has the advantage of being possible 
to implement without imbalance, and requires only local 
thermal equilibrium. Moreover fluctuations can be com- 
pared with susceptibilites determined from the equation 
of state to perform absolute, mo del- independent ther- 
mometry for strongly interacting systems [30]. 

In conclusion, we have demonstrated a new technique 
to determine spin susceptibilities of ultracold atomic 
gases using speckle imaging. We have validated and 
calibrated this technique using an ideal Fermi gas and 
applied it to a strongly interacting Fermi gas in the 
BEC-BCS crossover. This technique is directly appli- 
cable to studying pairing and magnetic ordering of two- 
component gases in optical lattices. 
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